Intelligent Modeling and Multi-Response Optimization of AWJC on Fiber Intermetallic Laminates through a Hybrid ANFIS-Salp Swarm Algorithm

The attainment of intricate part profiles for composite laminates for end-use applications is one of the tedious tasks carried out through conventional machining processes. Therefore, the present work emphasized hybrid intelligent modeling and multi-response optimization of abrasive waterjet cutting (AWJC) of a novel fiber intermetallic laminate (FIL) fabricated through carbon/aramid fiber, reinforced with varying wt% of reduced graphene oxide (r-GO) filled epoxy resin and Nitinol shape memory alloy as the skin material. The AWJC experiments were performed by varying the wt% of r-GO (0, 1, and 2%), traverse speed (400, 500, and 600 mm/min), waterjet pressure (200, 250, and 300 MPa), and stand-off distance (2, 3, and 4 mm) as the input parameters, whereas kerf taper (Kt) and surface roughness (Ra) were considered as the quality responses. A hybrid approach of a parametric optimized adaptive neuro-fuzzy inference system (ANFIS) was adopted through three different metaheuristic algorithms such as particle swarm optimization, moth flame optimization, and dragonfly optimization. The prediction efficiency of the ANFIS network has been found to be significantly improved through the moth flame optimization algorithms in terms of minimized prediction errors, such as mean absolute percentage error and root mean square error. Further, multi-response optimization has been performed for optimized ANFIS response models through the salp swarm optimization (SSO) algorithm to identify the optimal AWJC parameters. The optimal set of parameters, such as 1.004 wt% of r-GO, 600 mm/min of traverse speed, 214 MPa of waterjet pressure, and 4 mm of stand-off distance, were obtained for improved quality characteristics. Moreover, the confirmation experiment results show that an average prediction error of 3.38% for Kt and 3.77% for Ra, respectively, were obtained for SSO, which demonstrates the prediction capability of the proposed optimization algorithm.


Introduction
In recent times, composite materials are found to be an excellent alternative for structural, mechanical, and automotive applications instead of commercial materials such as metals and plastics because of the versatility in their properties, such as relative stiffness, high-strengthto-weight ratio, durability, rust, and corrosive resistance [1,2]. Fiber metallic laminates (FMLs) strongly impacting the chosen response characteristics is abrasive particle size. An experimental investigation was conducted by Hutyrová et al. [9] to investigate the surface topography of AWJ-machined wood plastic composites. According to their findings, AWJC can be utilized to effectively machine plastic composites without melting the matrix materials. For the assessment of the AWJC performance of jute/polyester composites, Kalirasu et al. [10] used an analytical method called response surface methodology and multiobjective optimization on the basis of ratio analysis (RSM-MOORA). They investigated how stand-off distance, cutting speed, and jet pressure affected the kerf taper and surface roughness of machined surfaces. The findings show that the proposed method can be used effectively for fiber-reinforced composites up to a maximum thickness of 3 mm. Using the response surface approach, K. Balamurugan et al. [11] investigated the AWJC characteristics of a composite made of lanthanum phosphate and yttria. They considered the impact of kerf taper, surface quality, and material removal rate as well as cutting speed, stand-off distance, and water pressure as input factors. Their findings showed that cutting speed has a negative impact on surface quality whereas water pressure has a positive impact on kerf taper and material removal rate. To examine the surface and kerf characteristics of AWJC of stacked Ti and carbon fibre reinforced plastics (CFRP) FMLs, Pahuja et al. [12] carried out experimental and statistical analyses. They asserted that the top skin is where micro buckling and fracture start to cause the erosion mechanism of FMLs. Additionally, the designs of the metal skin and polymer composite have a big impact on the FML interface failures. The effect of AWJC factors on the kerf quality of machined Ti/CFRP composite stacks has been analytically investigated by Ramulu [13]. They created an empirical model to anticipate the depth of penetration and material removal mechanism caused by waterjet pressure, and they discovered that the model can be successfully utilized to control AWJC parameters to produce cutting zones with no defects. On the qualitative traits of Ti/CFRP stacked composite laminates, a drilling operation through AWJ machining was carried out by Alberdi et al. [14]. Their findings showed that the layout of stacks, water pressure, and traverse speed of the focusing tube have a significant impact on the taper of drilling and surface quality.
Recently, few researchers have utilized intelligent modeling such as fuzzy logic, artificial neural network (ANN), adaptive neuro-fuzzy inference system (ANFIS), and metaheuristic optimization algorithms to obtain the optimal cutting conditions of the AWJC process [15][16][17][18][19]. However, none of the researchers have hybridized the metaheuristic algorithms to optimize ANFIS parameters to correlate the relationship between the input and response characteristics of the AWJC process, especially for fiber metallic laminates made of nickel titanium superalloy (Nitinol) as the skin material. A titanium dioxide-based protective surface layer protects nitinol from corrosion. Nitinol's biocompatibility and corrosion resistance are excellent when properly processed, and it provides the flexibility and strength required for many applications. When Nitinol is combined with carbon/aramid-based epoxy prepregs, it results in excellent structural stability with high impact strength.
The emphasis of the current work is on the AWJC of a novel fiber intermetallic laminate (FIL) made of sheets of Nitinol shape memory alloy and reduced graphene oxide (r-GO) filled epoxy prepregs embedded with carbon and aramid fibers, which are alternately interlaced. Among the available fabrication techniques, the vacuum resin infusion process has been utilized for fabricating the FILs due to their significant advantages such as high consistency, repeatability, minimized styrene emissions and a cleaner process than traditional FRP processing. The kerf taper and surface roughness were considered as independent response characteristics for the investigation. The AWJC experiments were carried out by varying processing parameters such as traverse speed, waterjet pressure, and stand-off distance as the independent variables with the different wt% of r-GO fillers in the FILs. Hybrid intelligent modeling of a parametric tuned ANFIS through metaheuristic algorithms was performed for correlating the relationship between the input parameters and output responses. Finally, the optimized ANFIS network was optimized through the salp swarm optimization (SSO) algorithm for obtaining optimal AWJC parameters.

Methodologies
The proposed work aims to adopt a hybrid approach of metaheuristic algorithm tuned ANFIS intelligent modeling of AWJC parameters and multi-response optimization through the salp swarm optimization algorithm. The flowchart shown in Figure 1 indicates the proposed methodology to achieve the modeling and optimization of the AWJC process. In the first phase, the ANFIS prediction model was proposed in order to correlate the AWJC parameters and output responses. The ANFIS model parameters were further optimized through metaheuristic algorithms such as particle swarm optimization (PSO), moth flame optimization (MFO), and dragonfly optimization (DFO), and the best algorithm was selected based on the obtained prediction errors such as root-mean-square error (RMSE) and mean absolute percentage error (MAPE). In the second stage, the ANFIS model was trained based on the optimal parameters and the prediction performance was evaluated with the actual response values. In the final stage, the optimized ANFIS network was considered as the objective function for multi-response optimization and the salp swarm optimization algorithm was conceived for obtaining the optimal AWJC parameters. The proposed methodologies such as response surface methodology (RSM), ANFIS, and the SSO algorithm were discussed in subsequent sections.

Response Surface Methodology
In order to develop empirical models without sacrificing their prediction ability, advanced manufacturing processes are adopting novel designs of experimental approaches to reduce the number of experimental trials. Since the AWJC process involves several uncertain processing conditions that significantly influence the quality features of the processed components, a competent experiment design approach, namely the Box-Behnken design (BBD) of response surface methodology, has been used for designing the viable number of experiments. In this present study, four parameters such as the addition of r-GO, traverse speed, waterjet pressure, and stand-off distance of each of the three levels have been considered for conducting the experiments with five replications for predicting the output responses such as kerf taper and surface roughness. The second order empirical models can be developed through the proposed BBD approach, which has been expressed as follows (Equation (1)): where  signifies statistical distribution error, i C , ij C , and ii C indicates the linear, interaction, and quadratic coefficients of the parameters, respectively, R indicates the response value and k indicates the number of dependent variables.

Adaptive Neuro-Fuzzy Inference System (ANFIS)
ANFIS is a hybrid machine learning intelligence model, which comprises the advantages of both fuzzy logic and neural network for precise mapping of the correlation that exists between the input and output parameters of complex real-world problems [20]. In ANFIS, the Takagi-Sugeno fuzzy inference system has been optimized by the highly interconnected neural network models, which enables the network to be more reliable for predicting the relationship between uncertain variables associated with any system. The general ANFIS network architecture consists of five major steps, as follows: input fuzzification, implication, normalization, defuzzification, and output layers. The detailed stepby-step explanation of how the ANFIS network has been effectively utilized for intelligent modeling of the advanced manufacturing process is explained by the previous works [21,22].

Response Surface Methodology
In order to develop empirical models without sacrificing their prediction ability, advanced manufacturing processes are adopting novel designs of experimental approaches to reduce the number of experimental trials. Since the AWJC process involves several uncertain processing conditions that significantly influence the quality features of the processed components, a competent experiment design approach, namely the Box-Behnken design (BBD) of response surface methodology, has been used for designing the viable number of experiments. In this present study, four parameters such as the addition of r-GO, traverse speed, waterjet pressure, and stand-off distance of each of the three levels have been considered for conducting the experiments with five replications for predicting the output responses such as kerf taper and surface roughness. The second order empirical models can be developed through the proposed BBD approach, which has been expressed as follows (Equation (1)): where ω signifies statistical distribution error, C i , C ij , and C ii indicates the linear, interaction, and quadratic coefficients of the parameters, respectively, R indicates the response value and k indicates the number of dependent variables.

Adaptive Neuro-Fuzzy Inference System (ANFIS)
ANFIS is a hybrid machine learning intelligence model, which comprises the advantages of both fuzzy logic and neural network for precise mapping of the correlation that exists between the input and output parameters of complex real-world problems [20]. In ANFIS, the Takagi-Sugeno fuzzy inference system has been optimized by the highly interconnected neural network models, which enables the network to be more reliable for predicting the relationship between uncertain variables associated with any system. The general ANFIS network architecture consists of five major steps, as follows: input fuzzification, implication, normalization, defuzzification, and output layers. The detailed step-by-step explanation of how the ANFIS network has been effectively utilized for intelligent modeling of the advanced manufacturing process is explained by the previous works [21,22].

Salp Swarm Optimization (SSO) Algorithm
Salp swarm optimization (SSO) is a novel swarm-based metaheuristic algorithm that mimics the navigating and foraging mechanism of salps in the oceans [23]. In the SSO algorithm, both leader and follower salps form the salp chain. The salp at the top of the salp chain is the leader, and the other salps are regarded as the followers. Half of the salps are chosen as leaders in order to increase the population diversity of the algorithm and the capacity to exit the local optima. The implementation flowchart of the proposed SSO algorithm for solving this engineering problem has been depicted in Figure 2. The step-by-step procedure of the SSO algorithm for obtaining optimal AWJC parameters is as follows: Step 1: Initialization of salp population A set of parameters P ij such as wt% of r-GO, TS, WP, and SOD are generated within lower bound and upper bound values as per their selected levels using the following relation (Equation (2)).
where lb j and ub j are the lower and upper bound value of j th parameter, and a rand-random value between 0 and 1.
Step 2: Evaluation of fitness of each Salp The fitness values of each salp, i.e., responses such as Kt and Ra, are calculated using the optimized ANFIS Model. The non-dominated salps, based on their dual objectives, are computed and the best salp, based on the crowding distance, is assumed as food source F 1j and stored in a file. The archive size is updated within its maximum size.
Step 3: Update the leader position and the other salps' position The value of the squared exponential covariance variable (C1), which defined the position of the leader salp, is calculated using the following relation (Equation (3)).
where it is the current iteration number and max_it is the maximum iteration number or stopping criteria. The leader salp's position is updated through the following relation (Equations (4) and (5)).
End where c 2 and c 3 are random values between 0 and 1, lb j and ub j are the lower and upper bound value of j th parameter, and F 1j is the position of food source. Further, the position of follower salps is described as follows (Equation (6)): Finally, the updated position of salps is verified, whether they are within their bounds, i.e., lb j and ub j .
Step 4: The initial position of salps is updated with the new updated position using the entire replacement strategy.
Step 5: Starting from step 2 to step 4, these steps are repeated until the stopping criteria.
Step 6: The dual objective solutions from the archive are converted into single objective solutions by implementing Deng's method. The order preference of the solution is made based on how high the value is in the overall performance index.
Step 7: The optimum process parameter (wt% of r-GO, TS, WP, and SOD) and its responses (Kt and Ra) are displayed based on the higher overall performance index value.
Materials 2022, 15, x FOR PEER REVIEW 6 of 20 Step 5: Starting from step 2 to step 4, these steps are repeated until the stopping criteria.
Step 6: The dual objective solutions from the archive are converted into single objective solutions by implementing Deng's method. The order preference of the solution is made based on how high the value is in the overall performance index.
Step 7: The optimum process parameter (wt% of r-GO, TS, WP, and SOD) and its responses (Kt and Ra) are displayed based on the higher overall performance index value.

Experimental Procedure and Measurements
The FILs are fabricated through a vacuum assisted resin infusion process by utilizing Nitinol shape memory alloy sheet as the skin material, bidirectional weaving patterned carbon and aramid fibers as the prepreg fiber materials, and reduced graphene oxide (r-GO), a derivative of graphene, which has been prepared through the modified hummer's method and has been supplemented as nano filler in the matrix material. The r-GO nano filled epoxy matrix is supplied through a resin inlet port of the infusion process in order to fabricate FILs with alternatively stacked woven carbon and aramid fibers with Nitinol shape memory alloy foils (NiTi/[carbon/epoxy/aramid]8/NiTi) as a skin material. The final thickness of fabricated FILs is upheld to approximately 3.5 mm.

Experimental Procedure and Measurements
The FILs are fabricated through a vacuum assisted resin infusion process by utilizing Nitinol shape memory alloy sheet as the skin material, bidirectional weaving patterned carbon and aramid fibers as the prepreg fiber materials, and reduced graphene oxide (r-GO), a derivative of graphene, which has been prepared through the modified hummer's method and has been supplemented as nano filler in the matrix material. The r-GO nano filled epoxy matrix is supplied through a resin inlet port of the infusion process in order to fabricate FILs with alternatively stacked woven carbon and aramid fibers with Nitinol shape memory alloy foils (NiTi/[carbon/epoxy/aramid]8/NiTi) as a skin material. The final thickness of fabricated FILs is upheld to approximately 3.5 mm.
The cutting experiments were performed on three different categories of FILs fabricated with 0, 1, and 2 wt% of r-GO nano fillers through a high-precision cantilever type computer numerical controlled AWJC system (S3015, Waterjet Germany, Chennai, India). The AWJ machine is equipped with a high-pressurized intensifier pump of operating pressure ranging up to 450 MPa and a carbide focusing tube of 0.76 mm diameter to achieve efficient machining operations. The machining parameters, such as the abrasive flow rate of 200g/min, the garnet abrasive particle size of~177µm (80 mesh), the orifice diameter of 0.276 mm and the jet impact angle of 90 • , are constantly maintained throughout the experiments. As per the RSM-BBD experimental design strategy for four dependent variables with three levels such as r-GO addition of 0, 1, and 2 wt%, traverse speed of 400, 500, and 600 mm/min, waterjet pressure of 200, 250, and 300 MPa, and stand-off distance of 2, 3, and 4 mm, a total of 29 straight parallel slots were made in the different FILs with a cutting length of 50 mm for the investigation of Kt and Ra. The range of AWJC parameters was selected by performing extensive preliminary studies. The typical AWJC experimental setup and processed FIL specimens is depicted in Figure 3. The AWJC system consists of abundant process-related dependent variables, whereas the most significant variables such as TS, WP, and SOD, suggested by previous researchers, were selected for the present investigation. The kerf taper at jet entry and exit of the machined FILs were computed using a METZ 1395 (Metzer Instruments, Mathura, India) tabletop microscope and the surface quality at the kerf side surface was measured using a non-contact three-dimensional profilometer (Talysurf CCI Lite, Taylor & Hobson, Leicester, UK). The computational errors were minimized by measuring Kt and Ra thrice and their mean values was considered as responses. The schematic representation of measuring Ra and Kt has been shown in Figure 4. The response features obtained from the AWJC of FILs are further considered for statistical validation and optimization, as presented in Table 1.
The cutting experiments were performed on three different categories of FILs fabricated with 0, 1, and 2 wt% of r-GO nano fillers through a high-precision cantilever type computer numerical controlled AWJC system (S3015, Waterjet Germany, Chennai, India). The AWJ machine is equipped with a high-pressurized intensifier pump of operating pressure ranging up to 450 MPa and a carbide focusing tube of 0.76 mm diameter to achieve efficient machining operations. The machining parameters, such as the abrasive flow rate of 200g/min, the garnet abrasive particle size of ~177μm (80 mesh), the orifice diameter of 0.276 mm and the jet impact angle of 90°, are constantly maintained throughout the experiments. As per the RSM-BBD experimental design strategy for four dependent variables with three levels such as r-GO addition of 0, 1, and 2 wt%, traverse speed of 400, 500, and 600 mm/min, waterjet pressure of 200, 250, and 300 MPa, and stand-off distance of 2, 3, and 4 mm, a total of 29 straight parallel slots were made in the different FILs with a cutting length of 50 mm for the investigation of Kt and Ra. The range of AWJC parameters was selected by performing extensive preliminary studies. The typical AWJC experimental setup and processed FIL specimens is depicted in Figure 3. The AWJC system consists of abundant process-related dependent variables, whereas the most significant variables such as TS, WP, and SOD, suggested by previous researchers, were selected for the present investigation. The kerf taper at jet entry and exit of the machined FILs were computed using a METZ 1395 (Metzer Instruments, Mathura, India) tabletop microscope and the surface quality at the kerf side surface was measured using a non-contact threedimensional profilometer (Talysurf CCI Lite, Taylor & Hobson, Leicester, UK). The computational errors were minimized by measuring Kt and Ra thrice and their mean values was considered as responses. The schematic representation of measuring Ra and Kt has been shown in Figure 4. The response features obtained from the AWJC of FILs are further considered for statistical validation and optimization, as presented in Table 1.   The cutting experiments were performed on three different categories of FILs fabricated with 0, 1, and 2 wt% of r-GO nano fillers through a high-precision cantilever type computer numerical controlled AWJC system (S3015, Waterjet Germany, Chennai, India). The AWJ machine is equipped with a high-pressurized intensifier pump of operating pressure ranging up to 450 MPa and a carbide focusing tube of 0.76 mm diameter to achieve efficient machining operations. The machining parameters, such as the abrasive flow rate of 200g/min, the garnet abrasive particle size of ~177μm (80 mesh), the orifice diameter of 0.276 mm and the jet impact angle of 90°, are constantly maintained throughout the experiments. As per the RSM-BBD experimental design strategy for four dependent variables with three levels such as r-GO addition of 0, 1, and 2 wt%, traverse speed of 400, 500, and 600 mm/min, waterjet pressure of 200, 250, and 300 MPa, and stand-off distance of 2, 3, and 4 mm, a total of 29 straight parallel slots were made in the different FILs with a cutting length of 50 mm for the investigation of Kt and Ra. The range of AWJC parameters was selected by performing extensive preliminary studies. The typical AWJC experimental setup and processed FIL specimens is depicted in Figure 3. The AWJC system consists of abundant process-related dependent variables, whereas the most significant variables such as TS, WP, and SOD, suggested by previous researchers, were selected for the present investigation. The kerf taper at jet entry and exit of the machined FILs were computed using a METZ 1395 (Metzer Instruments, Mathura, India) tabletop microscope and the surface quality at the kerf side surface was measured using a non-contact threedimensional profilometer (Talysurf CCI Lite, Taylor & Hobson, Leicester, UK). The computational errors were minimized by measuring Kt and Ra thrice and their mean values was considered as responses. The schematic representation of measuring Ra and Kt has been shown in Figure 4. The response features obtained from the AWJC of FILs are further considered for statistical validation and optimization, as presented in Table 1.

Statistical Analysis of Developed Mathematical Models
The adequacy and statistical significance of the developed mathematical models and performed experimental investigations were assessed through multi-analysis of variance (ANOVA) and Anderson-Darling normality tests. The experimental work considered the most significant AWJC parameters that influence the quality features of processed components. Therefore, the influence of individual parameters on the selected responses should be assessed through a systematic analysis for improving the cut quality characteristics. The ANOVA results for selected response features such as Kt and Ra have been mentioned in Tables 2 and 3. From the results, the coefficient of determination values for Kt and Ra is achieved at 99.1% and 99.2%, respectively. Moreover, the insignificant terms at a 95% confidence interval were eliminated from the regression modeling through the backward elimination procedure. The lack-of-fit values were also found to be significant within the selected range of processing parameters, which signifies the statistical significance of the conducted experiments. In addition, the statistical analysis (Anderson-Darling normality test) is performed to validate the obtained solutions at a 95% confidence level ( Figure 5). The results of statistical analysis also revealed that the obtained solutions are statistically significant and are normally scattered against actual values, which confirms the accuracy of conducted AWJC experiments. are statistically significant and are normally scattered against actual values, which confirms the accuracy of conducted AWJC experiments.

Intelligent Modeling of AWJC Quality Characteristics through Optimized ANFIS Network
Due to the uncertainty in non-traditional machining processes such as abrasive waterjet machining, the correlation that exists between the input parameters and the output characteristics is extremely non-linear. As a result, an efficient approach to modeling the relationship between the parameters is critical in order to improve the predictability

Intelligent Modeling of AWJC Quality Characteristics through Optimized ANFIS Network
Due to the uncertainty in non-traditional machining processes such as abrasive waterjet machining, the correlation that exists between the input parameters and the output characteristics is extremely non-linear. As a result, an efficient approach to modeling the relationship between the parameters is critical in order to improve the predictability of such processes. The correlation between the AWJC parameters such as traverse speed, waterjet pressure, and stand-off distance, as well as r-GO addition in FML on response features such as Kt and Ra was established in this work using the ANFIS intelligent modeling approach. To develop the proposed hybrid intelligent model, a customized code was written in the MATLAB 2021b™ environment. The proposed work utilized a sugeno-type sub-clustering "genfis2" ANFIS model for appropriate modeling of the inference system. Since the ANFIS network parameters, such as premise and consequent variables, are generally selected based on the choice of the user and/or a trial-and-error approach, this may lead to an inefficient model with a quantum of prediction errors. Therefore, an efficient tuning approach is necessitated to optimize these parameters to improve the prediction capability of the ANFIS network.
In the present work, three well-established metaheuristic algorithms namely PSO, MFO, and DFO have been efficiently utilized for tuning the ANFIS parameters in order to improve the learning capability by minimizing the prediction error. In general, the behavior of the subtractive clustering ANFIS model is influenced by the RADII, squelch factor, accept ratio, and rejection ratio. The RADII clustering center, which ranges from 0 to 1, represents the influencing ranges of each dependent and predictor variable. In several cases, it was discovered that a smaller cluster RADII is a better option for obtaining better prediction results. Therefore, a minimum range of RADII values between 0.13 and 0.5 were considered in the current work. The neighborhood cluster center was calculated by multiplying a chosen range of RADII values by the quash factor. Additionally, the volume of data required for effective training and performance testing of the FIS model is considered. The default values for the acceptance and rejection ratios are frequently 0.5 and 0.15, respectively. For each response, the RADII, squelch factor, number of training sessions, and data verification are additional factors that must be taken into consideration when building the best FIS model.
In the first step of ANFIS hybrid modeling, the control elements of the AWJC process, such as traverse speed, waterjet pressure, and stand-off distance with wt% of r-GO, are specified as input factors, whereas Kt and Ra are thought of as output factors. Furthermore, to increase the trained network's accuracy while reducing prediction errors, FIS parameters such as cluster radius, quash factor, and proportion of data used to train the network are selected. Table 4 contains a list of the ANFIS variable levels considered for the network's initial training with selected metaheuristic algorithms. The fuzzy rules are designed to cluster the chosen processing variables into a variety of values by combining two or more membership functions. Extensive membership functions must be generated in order to develop a rule-based correlation between the processing variables and response characteristics because there are numerous process parameters. As a result, this research utilizes the subtractive fuzzy clustering approach. Due to their smoothness and concise notation when estimating the responses, Gaussian-shaped membership functions (MFs) are preferred over a variety of other available membership functions (MFs). The key contributing parameters for PSO, MFO, and DFO for optimizing the ANFIS parameters have been mentioned in Tables 5 and 6. These parameters, which are the ideal training parameters for the current investigation, were obtained through multiple iterations of the trial-and-error method.

Parameters Value/Equation
No. of dragonflies (nd) The proposed algorithms were executed several times for improving the prediction capability of the ANFIS network by minimizing the MAPE and RMSE during the training of the network. The algorithms were simultaneously executed for 100 iterations in order to obtain the optimal ANFIS parameters. The convergence plots of DFO, MFO, and PSO algorithms as shown in Figure 6a-d, which indicates the correlation between the number of iterations and their consequences on the ANFIS outputs such as MAPE and RMSE for the response features such as Kt and Ra. The convergence plots show that the MFO algorithm has converged effectively and produced minimized MAPE and RMSE values when compared with the other two optimization algorithms. In addition, the statistical analysis results (normal distribution of data) have also proved that the MFO algorithm is efficient in the training of the proposed ANFIS network parameters at a 95% confidence interval, with an improved Anderson squared value of 0.72 for MAPE and 0.92 for RMSE, respectively (Figures 7 and 8

Achieve size 100
The proposed algorithms were executed several times for improving the prediction capability of the ANFIS network by minimizing the MAPE and RMSE during the training of the network. The algorithms were simultaneously executed for 100 iterations in order to obtain the optimal ANFIS parameters. The convergence plots of DFO, MFO, and PSO algorithms as shown in Figure 6a-d, which indicates the correlation between the number of iterations and their consequences on the ANFIS outputs such as MAPE and RMSE for the response features such as Kt and Ra. The convergence plots show that the MFO algorithm has converged effectively and produced minimized MAPE and RMSE values when compared with the other two optimization algorithms. In addition, the statistical analysis results (normal distribution of data) have also proved that the MFO algorithm is efficient in the training of the proposed ANFIS network parameters at a 95% confidence interval, with an improved Anderson squared value of 0.72 for MAPE and 0.92 for RMSE, respectively (Figures 7 and 8).       The optimal ANFIS parameters and their corresponding error values obtained for 30 epochs through the proposed optimization algorithms have been illustrated in Table 7. Among the available RMSE and MAPE values, MFO produced values are considered minimal for both the responses such as Kt and Ra. The optimized results obtained through MFO, 80% of the experimental data, i.e., 23 data, were considered for training the network, whereas the remaining six experimental data were taken to evaluate the prediction ability of the trained network. The objective of ANFIS intelligent modeling is to minimize the error between the training and checking data for each response. From the table, the training and checking errors obtained through the MFO algorithm for Kt is 0.0056% and 0.016%, and for Ra is 0.0088% and 0.0231%, respectively. The prediction error between the actual data and ANFIS predicted data has been graphically represented in Figure 9a,b. These findings suggest that the error values between the training and checking data are found to be minimal. As a result, this proposed trained ANFIS approach could be used to accurately predict AWJC parameters with fewer prediction errors. The optimal ANFIS parameters and their corresponding error values obtained for 30 epochs through the proposed optimization algorithms have been illustrated in Table 7. Among the available RMSE and MAPE values, MFO produced values are considered minimal for both the responses such as Kt and Ra. The optimized results obtained through MFO, 80% of the experimental data, i.e., 23 data, were considered for training the network, whereas the remaining six experimental data were taken to evaluate the prediction ability of the trained network. The objective of ANFIS intelligent modeling is to minimize the error between the training and checking data for each response. From the table, the training and checking errors obtained through the MFO algorithm for Kt is 0.0056% and 0.016%, and for Ra is 0.0088% and 0.0231%, respectively. The prediction error between the actual data and ANFIS predicted data has been graphically represented in Figure 9a,b. These findings suggest that the error values between the training and checking data are found to be minimal. As a result, this proposed trained ANFIS approach could be used to accurately predict AWJC parameters with fewer prediction errors.

Influence of Process Parameters on Quality Characteristics
The kinetic energy of the waterjet and its associated variables, as well as the thickness of the substrate materials, have a significant impact on the cutting quality in the AWJC process. The primary cause of an irregular kerf taper ratio, which lowers the calibers of machined components, is the energy depletion of the waterjet along the thickness of the substrate. To achieve the desired cutting quality with minimized taper, appropriate AWJC variables must be chosen and their effects on kerf taper must be investigated. Moreover, to determine how the machined components will interact with their adjacent components and various working environments, it is crucial to investigate the surface quality of abrasive waterjet machined materials, especially layered composites.
The influence of various AWJC parameters and the inclusion of r-GO with varying weight percentages on the quality characteristics of processed FILs at varying cutting conditions are described in the ANFIS surface plots (Figure 10a-d). The influence of waterjet pressure (WP) and wt% of r-GO on Kt is represented in the 3D surface plot (Figure 10a). The plot reveals that Kt linearly increases with an increase in WP from 200 bar to 300 bar, whereas a slight variation in Kt has been found when the addition of r-GO has been increased from 0 to 2 wt%. When the waterjet pressure is increased, the FILs may be impinged upon by the kinetic energy and flow turbulence of the water-jet stream, creating a wider slot. In addition, interparticle collision is increased by flow turbulence and waterjet expansion, which widens the kerf surface's taper [12,24]. The consequences of traverse speed (TS) and stand-off distance (SOD) on Kt have been demonstrated in Figure 10b. The plot shows that Kt is linearly widening with an increase in SOD from 2 mm to 4 mm and TS from 400 mm/min to 600 mm/min. This might be because of a narrowing of the effective waterjet caused by higher SOD, which causes the waterjet to flow away from the focusing tube and contract into abrasive particles. As a result, the cutting and penetrating zones of abrasive particles in a waterjet are reduced in energy. The kerf on the bottom cut surface is consequently smaller than the kerf on the top cut surface. Kt was therefore found to be greater at higher SOD levels [25].

Influence of Process Parameters on Quality Characteristics
The kinetic energy of the waterjet and its associated variables, as well as the thickness of the substrate materials, have a significant impact on the cutting quality in the AWJC process. The primary cause of an irregular kerf taper ratio, which lowers the calibers of machined components, is the energy depletion of the waterjet along the thickness of the substrate. To achieve the desired cutting quality with minimized taper, appropriate AWJC variables must be chosen and their effects on kerf taper must be investigated. Moreover, to determine how the machined components will interact with their adjacent components and various working environments, it is crucial to investigate the surface quality of abrasive waterjet machined materials, especially layered composites.
The influence of various AWJC parameters and the inclusion of r-GO with varying weight percentages on the quality characteristics of processed FILs at varying cutting conditions are described in the ANFIS surface plots (Figure 10a-d). The influence of waterjet pressure (WP) and wt% of r-GO on Kt is represented in the 3D surface plot (Figure 10a). The plot reveals that Kt linearly increases with an increase in WP from 200 bar to 300 bar, whereas a slight variation in Kt has been found when the addition of r-GO has been increased from 0 to 2 wt%. When the waterjet pressure is increased, the FILs may be impinged upon by the kinetic energy and flow turbulence of the water-jet stream, creating a wider slot. In addition, interparticle collision is increased by flow turbulence and water-jet expansion, which widens the kerf surface's taper [12,24]. The consequences of traverse speed (TS) and stand-off distance (SOD) on Kt have been demonstrated in Figure 10b. The plot shows that Kt is linearly widening with an increase in SOD from 2 mm to 4 mm and TS from 400 mm/min to 600 mm/min. This might be because of a narrowing of the effective waterjet caused by higher SOD, which causes the waterjet to flow away from the focusing tube and contract into abrasive particles. As a result, the cutting and penetrating zones of abrasive particles in a waterjet are reduced in energy. The kerf on the bottom cut surface is consequently smaller than the kerf on the top cut surface. Kt was therefore found to be greater at higher SOD levels [25].
The combinatory effects of the AWJC process parameters on Ra have been represented by the ANFIS three-dimensional plots (Figure 10c,d). The influence of waterjet pressure and r-GO addition on surface quality has been presented in Figure 10c. It is perceived from the plot that Ra has been increased linearly with an increase in waterjet pressure, whereas the addition of r-GO nano fillers is found to be insignificant on surface quality. As the waterjet pressure increases, the energy fluctuation of the waterjet also intensifies. The high-pressure waterjet's erosion of abrasive particles will make the cut surface more uneven, which results in a rough-cut surface [26]. A three-dimensional surface plot (Figure 10d) shows the effect of the stand-off distance and traverse speed on the Ra of processed FIL surfaces. As can be seen from the plot, the Ra is significantly increasing, with an increase in both the processing variables from their low to high levels. When the stand-off distance rises, the waterjet expands before impacting the substrate, exposing the substrate to more extrinsic drag from the surrounding environment. As a result, the kinetic energy of the jet is reduced, causing the jet diameter to expand, resulting in a rough-cut surface [25,27].
sure, whereas the addition of r-GO nano fillers is found to be insignificant on surface quality. As the waterjet pressure increases, the energy fluctuation of the waterjet also intensifies. The high-pressure waterjet's erosion of abrasive particles will make the cut surface more uneven, which results in a rough-cut surface [26]. A three-dimensional surface plot (Figure 10d) shows the effect of the stand-off distance and traverse speed on the Ra of processed FIL surfaces. As can be seen from the plot, the Ra is significantly increasing, with an increase in both the processing variables from their low to high levels. When the stand-off distance rises, the waterjet expands before impacting the substrate, exposing the substrate to more extrinsic drag from the surrounding environment. As a result, the kinetic energy of the jet is reduced, causing the jet diameter to expand, resulting in a roughcut surface [25,27].

Multi-Response Optimization through Salp Swarm Algorithm
In this study, a hybrid approach of a parametric tuned ANFIS-salp swarm optimization algorithm has been used for the modeling and optimization of AWJC parameters during the processing of novel fiber intermetallic laminates. Based on the results of the proposed hybrid approaches, the following inferences were attained: The present investigation aims to optimize the AWJC parameters in order to improve the cut quality characteristics such as kerf taper and surface quality. In general, the multiresponse optimizations were performed in two ways: (i) the objective functions were combined into single objectives by assigning weight for each function to minimize the optimization complexity, and (ii) to obtain optimal results through pareto-optimal solutions of the deciding parameters [28]. In this study, the optimization was performed through a

Multi-Response Optimization through Salp Swarm Algorithm
In this study, a hybrid approach of a parametric tuned ANFIS-salp swarm optimization algorithm has been used for the modeling and optimization of AWJC parameters during the processing of novel fiber intermetallic laminates. Based on the results of the proposed hybrid approaches, the following inferences were attained: The present investigation aims to optimize the AWJC parameters in order to improve the cut quality characteristics such as kerf taper and surface quality. In general, the multi-response optimizations were performed in two ways: (i) the objective functions were combined into single objectives by assigning weight for each function to minimize the optimization complexity, and (ii) to obtain optimal results through pareto-optimal solutions of the deciding parameters [28]. In this study, the optimization was performed through a novel metaheuristic-based salp swarm optimization algorithm in order to obtain the non-dominated optimal solutions. The trained ANFIS outputs for each response have been considered as the objective functions for the multi-response optimization. The SSO algorithm was executed with the correlation of terms as indicated in Table 8. The implementation of SSO is simple and only needs a few parameters. Therefore, the algorithm was able to produce the best results possible with these minimal parameter settings. The optimization was executed by using the SSO algorithm with the parameter settings of 30 number of salps with a maximum of 100 iterations. During the optimization, the parameters of SSO were self-tuned for obtaining global optimal solutions that satisfy the goal of objective functions. The algorithms were executed twenty-nine times to attain ingenious pareto-optimal settings. From each execution, a set of pareto-optimal solutions were attained and the best one among the optimal front was selected through Deng's statistical approach [29]. The optimization performance of the SSO algorithm has been shown in the sample pareto-optimal front ( Figure 11). Among the obtained optimal parameter settings through each run, the best one was attained by implementing the Deng's statistical ranking approach. The solutions were ranked based on the Deng's value, as mentioned in Table 9. From the table, a higher Deng's value of 0.51507 was obtained at execution number 15 (As highlighted in bold) and their corresponding parameters, such as 1.004 wt% of r-GO, 600 mm/min of TS, 214 MPa of WP, and 4 mm of SOD, were global optimal parameter settings that provide an improved response characteristic value such as Kt of 2.067 • and Ra of 3.18 µm, respectively. In addition, a statistical analysis (Anderson-Darling normality test) was performed to validate the obtained solutions at a 95% confidence level ( Figure 12). The results of statistical analysis revealed that the obtained solutions are statistically significant and are normally scattered against actual values, which confirms the efficacy of the proposed optimization algorithm. novel metaheuristic-based salp swarm optimization algorithm in order to obtain the nondominated optimal solutions. The trained ANFIS outputs for each response have been considered as the objective functions for the multi-response optimization. The SSO algorithm was executed with the correlation of terms as indicated in Table 8. The implementation of SSO is simple and only needs a few parameters. Therefore, the algorithm was able to produce the best results possible with these minimal parameter settings. The optimization was executed by using the SSO algorithm with the parameter settings of 30 number of salps with a maximum of 100 iterations. During the optimization, the parameters of SSO were self-tuned for obtaining global optimal solutions that satisfy the goal of objective functions. The algorithms were executed twenty-nine times to attain ingenious pareto-optimal settings. From each execution, a set of pareto-optimal solutions were attained and the best one among the optimal front was selected through Deng's statistical approach [29]. The optimization performance of the SSO algorithm has been shown in the sample pareto-optimal front ( Figure 11). Among the obtained optimal parameter settings through each run, the best one was attained by implementing the Deng's statistical ranking approach. The solutions were ranked based on the Deng's value, as mentioned in Table 9. From the table, a higher Deng's value of 0.51507 was obtained at execution number 15 (As highlighted in bold) and their corresponding parameters, such as 1.004 wt% of r-GO, 600 mm/min of TS, 214 MPa of WP, and 4 mm of SOD, were global optimal parameter settings that provide an improved response characteristic value such as Kt of 2.067° and Ra of 3.18 μm, respectively. In addition, a statistical analysis (Anderson-Darling normality test) was performed to validate the obtained solutions at a 95% confidence level ( Figure 12). The results of statistical analysis revealed that the obtained solutions are statistically significant and are normally scattered against actual values, which confirms the efficacy of the proposed optimization algorithm. Figure 11. Sample pareto optimal front of SSO algorithm for optimized AWJC parameters. Figure 11. Sample pareto optimal front of SSO algorithm for optimized AWJC parameters.   Confirmation experiments were conducted to assess the rationality of the proposed optimization approach. The confirmation experiments were repeated thrice with the obtained optimal AWJC parameters through the SSO algorithm and their average values are presented in Table 10. The table shows that the correlation between predicted and experimentally measured response values is relatively good, with an average error of 3.38% for Confirmation experiments were conducted to assess the rationality of the proposed optimization approach. The confirmation experiments were repeated thrice with the obtained optimal AWJC parameters through the SSO algorithm and their average values are presented in Table 10. The table shows that the correlation between predicted and experimentally measured response values is relatively good, with an average error of 3.38% for Kt, and 3.77% for Ra, respectively. Therefore, the proposed SSO algorithm can be adequate for obtaining optimal processing conditions during AWJC of r-GO reinforced FILs for better cut quality characteristics.

Conclusions
In this study, a hybrid approach of a parametric tuned ANFIS-salp swarm optimization algorithm has been used for the modeling and optimization of AWJC parameters during the processing of novel fiber intermetallic laminates. Based on the results of the proposed hybrid approaches, the following inferences were attained:

•
The Box-Behnken design-based response surface methodology approach has been effectively utilized for designing the experimental trails with a viable number of experiments to minimize the time and cost. • The MFO algorithm was found to be outperformed in the training of the ANFIS network parameters with an average MAPE of 0.16838 for kerf taper and 0.10713 for surface roughness, and RMSE of 0.00895 for kerf taper and 0.01315 for surface roughness, respectively.

•
The hybrid MFO-ANFIS approach was found to be an efficient approach for intelligent modeling of the AWJC process to correlate its process parameters and response characteristics, with average training and testing errors of 0.0056% and 0.0163% for kerf taper, and 0.0088% and 0.0231% for surface roughness, respectively.

•
The optimal AWJC conditions for improved quality characteristics have been achieved through a metaheuristic salp swarm optimization algorithm, which are: 1.004 wt% of r-GO, 600 mm/min of TS, 214 MPa of WP, and 4 mm of SOD for an improved response characteristic value such as Kt of 2.067 • and Ra of 3.18 µm, respectively.

•
The results of confirmation studies revealed that the proposed optimization approach can be efficient in the prediction of optimal AWJC parameters with an error between the experimental and SSO algorithm predicted values of 3.38% for Kt and 3.77% for Ra, respectively.

Conflicts of Interest:
The authors declare no conflict of interest.